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Abstract 


The effects of the first nonlinear corrections to the DGLAP evolntion eqnations are stndied 
by using the recent HERA data for the structure function F 2 {x, Q‘^) of the free proton and the 
parton distributions from CTEQ5L and CTEQ6L as a baseline. By requiring a good fit to the 
HI data, we determine initial parton distributions at Qq = 1.4 GeV^ for the nonlinear scale 
evolution. We show that the nonlinear corrections improve the agreement with the F 2 {x,Q‘^) 
data in the region of x ~ 3 ■ 10“® and ~ 1.5 GeV^ without paying the price of obtaining a 
worse agreement at larger values of x and Q^. For the gluon distribution the nonlinear effects 
are found to play an increasingly important role at x < 10“^ and <10 GeV^, but rapidly 
vanish at larger values of x and Q^. Gonsequently, contrary to GTEQ6L, the obtained gluon 
distribution at = 1.4 GeV^ shows a power-like growth at small x. Relative to the GTEQ6L 
gluons, an enhancement up to a factor ~ 6 at x = 10“®, Ql = 1.4 GeV^ reduces to a negligible 
difference at Q^>10 GeV^. 
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1 Introduction 


Parton distribution functions (PDFs), fi{x,Q^), are needed for the computation of inclusive 
cross sections of hard, collinearly factorizable, processes in hadronic collisions. At sufficiently 
large values of the interaction scale and the momentum fraction x, where power corrections 
are negligible, the scale evolution of the PDFs is predicted quite accurately by the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations [1] derived from perturbative 
QGD. 

In the global analyses of the PDFs of the free proton, such as in Refs. [2, 3], the lowest 
order (LO), next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) parton 
distributions are extracted within the DGLAP framework using constraints from the measured 
cross sections of various hard processes and from the sum rules. In the proton case the procedure 
is well established: the initial distributions given at some initial scale Qq are hrst evolved to 
larger by using the DGLAP equations, then a comparison with the data is made over a wide 
range of x and Q^, after which the initial distributions are iterated until a good global £t to 
the data is obtained. The initial distributions are thus the non-perturbative input needed, the 
element that perturbative QGD cannot predict. The data from deeply inelastic lepton-proton 
scattering (DIS) play a key role in these analyses especially in the region of the smallest values 
of x, where the DIS data from ep collisions at DESY-HERA give the only constraints available. 
Recently, the HI collaboration at HERA [4] has measured the structure function F 2 {x,Q‘^) of 
the proton down to a; ~ 3 ■ 10“^ but still in the perturbatively accessible region > 1.5 GeV^. 
These data have been included in the recent global analyses by the MRST [2] and GTEQ [3] 
collaborations. 

In spite of the impressive success of the DGLAP approach, certain problems appear in the 
attempts to make the global fits to the HI data [4] as good as possible simultaneously in the 
region of > 4 GeV^ (“large” and in the region of 1.5 GeV^ < < 4 GeV^ (“small” Q^). 

In the recent NLO analysis MRST2001 [2] both regions are included and a good overall fit is 
found but with the expense of allowing for a negative gluon distribution. Although a negative 
contribution in the NLO gluon distribution is acceptable as long as the NLO cross sections 
remain positive, the interpretation of the PDFs as probability or number density distributions 
becomes obscured On the other hand, the GTEQ collaboration emphasizes the large-Q^ 
region in their global fits: e.g. in the sets GTEQ5 [6] and GTEQ6 [3] only the region > 4 
GeV^ of the DIS data is included in the fit and a very good agreement with the data is found. 
The agreement of the extrapolation to the small-Q^ region, however, becomes then worse. 
One is facing the problem of negative gluon distributions also in the NLO set GTEQ6M [3]: 
xg{x, is set to zero at the smallest values of x at scales < 1.69 GeV^. In LO, the negative 
gluon distributions, however, are not allowed. The quality of the LO fits to the HI data (see 
Table 1 in Sec. 3) reflects the problem of a simultaneous fit to the small- and large-Q^ regions: 
MRST2001 (GTEQ6) fits the small-Q^ (large-Q^) region better. 

The problems discussed above are very interesting as they can be a sign of a new QGD 
phenomenon: towards smaller values of x and (or) (but still 3> Aq^^j-,), gluon recombina¬ 
tion effects are expected to play an increasingly important role. These effects induce nonlinear 
power corrections to the DGLAP equations. First of the nonlinear terms have been calculated 

®See also [5] for a recent discussion of the PDFs not being probabilities. 
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by Gribov, Levin and Ryskin in [7] and by Mueller and Qiu in [8]. We shall refer to these 
corrections as the GLRMQ terms. 

Previous studies of the GLRMQ terms in the context of extracting the PDFs of the free 
proton can be found e.g. in [9]. Also other nonlinear evolution equations relevant at high 
gluon densities have been derived in the recent years [10], and the structure functions from DIS 
have been analysed in the context of saturation models [11]. In the present work, however, we 
shall adopt the framework of collinear factorization with universal PDFs, and search for the 
nonlinear GLRMQ corrections on top of the full DGLAP equations. This allows for a direct 
comparison of our results with those of the global DGLAP hts [2, 3]. In this manner we can 
also show more explicitly the need for nonlinear terms in the evolution equations. 

In the DGLAP evolution, the dependence of the sea quarks at small values of x is dictated 
by the gluon distribution; the larger xgi^XjQ"^), the faster the evolution of F 2 {x,Q'^), since 
in LO dF 2 {x,Q‘^)/d\ogQ‘^ ^ {10as/277r)xg{2x,Q'^) [12]. The GLRMQ terms slow down the 

evolution of gluons and sea quarks from the standard DGLAP behaviour (assuming the 
same starting distributions). Gonsequently, the logQ^ slopes of F 2 {x,Q‘^) of the HI data can 
be reproduced with a larger gluon distribution than that in the conventional DGLAP case. In 
this paper, we investigate this interdependence of the initial conditions at Ql = lA GeV^ and 
the effects of the GLRMQ terms by using the HI data as a baseline. We demonstrate that 
inclusion of the GLRMQ terms on top of the LO DGLAP evolution improves the agreement 
with the data in the small-a; and small-Q^ region while simultaneously maintains the good £t 
of the LO sets of GTEQ5 and GTEQ6 to the data at larger values of x and Q^. We show 
that the obtained small-a; gluon distribution can still have a power-like growth at the scale 

= 1.4 GeV^, leading to an enhancement of a factor ~ 6 relative to the GTEQ6L gluons at 
X = 10“®. We also show explicitly how the large deviations from GTEQ6L reduce to very small 
differences at > 10 GeV^. The size of the nonlinear terms, uncertainties and applicability 
region of the DGLAP-l-GLRMQ approach are discussed as well. 


2 Nonlinear evolution equations 


The GLRMQ corrections [7, 8] arise from fusion of two gluon ladders, and they modify the 
evolution equations of gluons as 


dxg{x,Q^) _ dxg{x,Q^) _ A ^ 2 ^( 2 )/ /.n 

aiogQ2 aiogQ2 DGLAP 2Q2^ yV U 

where the first term is the standard DGLAP result [1], linear in the PDFs. The 2-gluon density 
in the second term we model as 


x^G^‘^\x, Q'^) 


1 


[xg{x,Q^)f, 


( 2 ) 


with i? = 1 fm as the radius of the proton. The equation (1) thus becomes nonlinear in xg. 
In the scale evolution of the sea quark distributions the leading nonlinear correction from the 
gluon fusion appears as [8] 
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where the hrst term is again from the DGLAP equations and linear in the PDFs. Note that 
here we work under an approximation of neglecting contribution from the “high twist” gluon 
distribution G-in^x^Q"^) [8]. More discussion on this will be given later [13]. The evolution of 
the valence quark distributions remains unmodihed. 


3 Analysis and Results 


3.1 Effects of the GLRMQ corrections 

The emphasis of the CTEQ analysis is in the large-Q^ region where the nonlinearities should 
remain small. Since we wish to search for the nonlinear effects in the small-Q^ (small-x) region 
but also recover the DGLAP evolution at large the GTEQ distributions are an ideal baseline 
for our analysis. 

In Fig. 1, together with the HI data [4], we plot the dependence of the LO structure 
function F 2 {x,Q‘^) = J2q€^[xq{x,Q‘^) + xq{x,Q‘^)] computed from GTEQ5L [6, 14] for fixed 
values of x (dashed lines). The agreement with the HI data is clearly getting worse towards 
smaller values of x and Q^. This trend is clearly visible also in Fig. 2 (dashed line), where we 
quantify the quality of the ht in terms of 




n(xj) rpth/ _ pexp/ 

_ J] 


(4) 


divided by the cumulative number of data points, N{xk) = J 2 x^=o .2 a function of the x 

of the data. Here n{xj) refers to the number of the data points with the same value of a; = Xj. 
The x^/A^ computed at different regions of is in turn summarized in Table 1. Relative to 
the data [4] (which was not available at the time of GTEQ5) the computed logQ^ slopes of F 2 
are too large at small values of x and Q^. This is caused by a too large gluon component in 
the DGLAP evolution at small values of x. 

For showing how the GLRMQ terms slow down the scale evolution, we take the PDFs from 
the GTEQ5L parametrization at = 5 GeV^, and evolve these both downwards and upwards 
by using the nonlinear evolution equations. The results are plotted in Fig. 1 (solid lines). As 
seen in the hgure, at larger values of x and the nonlinear effects remain small and do not 
make the agreement with the data worse. In addition, as shown both by Fig. 1 and Fig. 2, 
at 5 ■ 10“®<a;< 10“^ the agreement with the data becomes quite good (i.e. x^/A^ remains 
constant), clearly improving the situation from the GTEQ5L case. 

Another interesting observation from Fig. 1 is that, with the trial initial conditions above, 
the nonlinear corrections at a; < 5 ■ 10“® and ~ 1 GeV^ obviously become too large and 
they start to dominate the computed evolution, causing negative logQ^ slopes for F 2 {x,Q^). 
Glearly, these are not supported by the data. In the computation we would then have entered 
a gluon saturation region where also other correction terms in the evolution equations should 
be included. This region is thus beyond the scope of the present DGLAP+GLRMQ approach. 
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Q" < 4.0 GeV" 
N = 29 

g" > 4.0 GeV" 
N = 104 

all g" 

N = 133 

GTEQ5L 

31.8 

1.18 

7.86 

GTEQ6L 

2.72 

0.93 

1.32 

MRST2001 

0.59 

2.06 

1.74 

This work: 

Set 1: Qc < GeV 

1.75 

0.96 

1.13 

Set 2a: Qc = L3 GeV 

1.58 

1.05 

1.17 

Set 2b: Qe = yiA GeV 

0.95 

0.86 

0.88 


Table 1: The /N of the fit of the computed LO F 2 to the HI data [4] in the small-Q^ and large-Q^ 
regions. Also the number of data points in each case is mentioned. 


The too strong nonlinearities can again be traced back to a too large gluon component at small 
X in the initial condition chosen. A more realistic gluon distribution which evolves according 
to the DGLAP+GLRMQ equations should obviously be smaller than GTEQ5L at small values 
of a; at = 5 GeV^. 

The HI data [4] is partly taken into account in the latest set GTEQ6L [3], and the situation 
is clearly improved from GTEQ5L. This is shown in Figs. 3 and 2 by the dotted-dashed curves. 
In Fig. 4 we have plotted the PDFs at = 1.4 GeV^ for several cases. Now, as shown by Fig. 
4, the gluon distributions (upper left panel) of GTEQ6L (dotted-dashed) at small x are smaller 
than those of GTEQ5L (dashed), causing more modest logQ^ slopes for F 2 at small x. Fig. 2 
indicates that the agreement of GTEQ6L with the HI data is excellent at a; > 10“^ and stays 
very good also at smaller x. Table 1 again expresses y^/A^ at the large-Q^ region (included in 
the GTEQ6 global £t) and in the small-Q^ region (not included in the GTEQ6 fit). Although a 
good agreement with the data is found, a closer look at Fig. 3 shows that again the agreement 
of the pure DGLAP result is getting worse towards the smallest values of x. 

We notice that, contrary to GTEQ5L, the GTEQ6L result in Fig. 3 lies above the HI data 
at the smallest values of x. If we now simply took the initial conditions from GTEQ6L at, 
say, = b GeV^, and evolved the distributions downwards by using the nonlinear equations 
(which would slow down the evolution), we would make the agreement with the data in the 
small-x region worse than the GTEQ6L result. The dotted curve in Fig. 2 shows the y^/A^ 
for such a run, notice the growing trend towards the smallest x. Encouraged, however, by the 
observations with the GTEQ5L-I-GLRMQ above, we wish to see whether we could find initial 
conditions at some Qq that would lead to at least the same or possibly even better agreement 
with the HI data as the DGLAP result from GTEQ6L. 
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3.2 New initial distributions 
3.2.1 Set 1 


We construct the initial distributions at an initial scale Q‘l = lA GeV^ by using the CTEQ5L 
and CTEQ6L sets as our guide. Following CTEQ5L we use = 192 MeV together with 

the one-loop expression of the strong coupling constant. The change of the value of Aqcd at 
the heavy-quark mass thresholds is taken into account. Throughout the study we use Qh = 
4.75 GeV for the 6-threshold, and the c-threshold will be discussed below. As seen in Fig. 2 
(sparsely dotted curve), a good agreement with the HI data is found at a; > 0.01 by taking the 
GTEQ6L PDFs at 5.0 GeV^ and evolving them down to 1.4 GeV^ according to the nonlinear 
equations. The GTEQ5L distributions evolved down to 1.4 GeV^ from 10.0 GeV^ and 3.0 
GeV^ with the GLRMQ corrections included, give a reasonably good agreement with the HI 
data at 10“"^ < 0.01 and at 10“^ < 10“^, correspondingly (not shown). A working initial 

condition can then be found by interpolating between these three results. Finally, we make a 
power-law £t to the interpolated gluon and light sea-quark distributions at small values of x. 
The x-slope of the small-x gluons is tuned to reproduce the measured logQ^ slopes of F 2 . This 
leads to an initial gluon distribution xg{x,Ql) = 3.64 ■ (0.01/a;)°'^® at small x. 

The initial conditions constructed in this way at = 1.4 GeV^ are shown in Fig. 4, labelled 
as our “set 1”. The obtained distributions are compared with the GTEQ5L and GTEQ6L PDFs 
at the same scale. Notice especially that our gluon distribution at a; = 10“® at this scale is a 
factor 6.4 larger than the one in GTEQ6L. In turn, the light sea-quark distributions at small x 
he in between GTEQ6L and GTEQ5L. Due to the procedure we have chosen, the valence-quark 
distributions (which evolve according to DGLAP) are practically the same as in GTEQ5L at 
a; <0.01 and GTEQ6L at a; >0.01. 

The DGLAP-I-GLRMQ evolution to higher scales then gives the result shown in Fig. 3 by 
the solid lines. The corresponding values of are again shown as a function of x in Fig. 

2 and its division into small-Q^ and large-Q^ regions in Table 1. We observe that while the 
agreement with the HI data at large values of x and is maintained practically as good as in 
GTEQ6L, the fit in the smallest-a;, small-Q^ region is indeed improving. This, together with 
the parton distributions at = 1.4 GeV^ for the nonlinear evolution, is the main result of this 
paper. At this point we should also emphasize that we have not attempted to make a global 
statistical analysis such as the one by GTEQ in order to further minimize the y^. After such a 
procedure, further improvement on the y^ could well be anticipated. 

Fig. 5 shows the dependence of the gluon distributions at fixed values of x as obtained 
in the DGLAP analyses GTEQ5L (dashed lines) and GTEQ6L (dotted-dashed) as well as from 
the GLRMQ-I-DGLAP evolution of the gluons of our set 1 (solid). Notice on one hand the large 
difference in the GTEQ5L and GTEQ6L sets, and on the other hand the slower evolution in 
the nonlinear case at small values of x and Q^. We also draw attention to the fact that in spite 
of the large (factor 6.4) difference at Ql = 1.4 GeV^ at a; = 10“^, our gluon distributions and 
the GTEQ6L gluons are in fact quite similar at > 5 GeV^. Goming back to the discussion 
related to Fig 1, we also observe from Fig. 5 that at = 5 GeV^ the gluons from our set 1 
are below the GTEQ5L gluons as anticipated. 
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3.2.2 Sets 2a and 2b 


A detail to study next is the c-quark contribution. The initial distributions in our set 1 were 
obtained on the basis of distributions which were evolved downwards from 3, 5 and 10 GeV^ 
according to the nonlinear evolution equations. The c-quarks we treat as massless quarks which 
experience the GLRMQ corrections as well. Therefore, the downwards evolution of the c-quarks 
is also slower than in the DGLAP case, and their distributions vanish only somewhat below 
the threshold scale Qc = 1.3 GeV of the sets GTEQ6L and GTEQ5L. This explains the fairly 
large c-distribution at Ql in our set 1. 

We construct two slightly different sets of PDFs with the same initial conditions at Qq = 
1.4 GeV^, where xc{x,Ql) = 0. For the first case, called “set 2a” we follow GTEQ6L and 
take Qc = 1.3 GeV as the c-threshold scale. For the other case, called “set 2b” we choose, 
Qc = Qo = Vl.4 GeV. In these initial conditions, the gluons are not modihed from set 1. 
However, in order to compensate for the loss of the initial c-quarks and to recover the agreement 
with the HI data, we enhance the other sea-quarks in the small-a; region simply by tuning their 
power in the power-law fit. The initial distributions for our sets 2 are shown in Fig. 4 (double 
dashed). After the nonlinear evolution to higher scales, the agreement of the set 2a with the 
HI data is practically as good as with our set 1 and GTEQ6L, as is shown in Figs. 3 and 2 
(double dashed lines). With the set 2b, the agreement becomes even better, as can be seen 
from Figs. 2 and 3 (short dashed lines). The same conclusion is suggested also by the x^/iV 
computed in the different regions of Q^ in Table 1. The GLRMQ corrections thus improve the 
hts to the data in the region of small x and Q^ without loosing the good hts at larger x and Q^. 
Regarding the sensitivity of the results to the c-threshold, we note that the differences between 
our sets 2a and 2b are quite small and could most probably be obtained also by keeping the 
c-threshold constant while tuning the gluon distributions at 0.001 < a; <0.01. This hne-tuning 
is, however, beyond the goal of this paper. 


4 Discussion and Conclusions 

It is important to keep in mind the uncertainties and limitations of the present approach. A 
constraint not addressed above but used in the global analysis of the PDFs, is momentum con¬ 
servation. Due to the nonlinear terms in Eqs. (1) and (3), some momentum from small values 
of X is lost in the evolution towards higher scales. Due to this, our sets 1 and 2 overestimate 
the total momentum at Q^ = 1.4 GeV^ by less than 2 %. By Q^ = 100 GeV^, some 2.6 % of 
the initial total momentum is lost. As the emphasis of our study is in the small-x behaviour 
of the PDFs, and as the violation remains small, we have not made an attempt to correct the 
obtained distributions for the momentum sum rule. 

In modelling the 2-gluon density in Eq. (1), we have set the effective radius parameter of 
the free proton to R = 1 fm. Depending on the transverse matter density prohle assumed for 
the free proton, some ~ 20 % uncertainty in R can be expected. The nonlinearities decrease 
with increasing R, so due to the interplay between the initial conditions and the scale evolution 
demonstrated above, a larger R would lead to a smaller enhancement of the small-x gluons. 
In order to properly estimate the uncertainty in the initial gluon distribution caused by the 
uncertainty in R, the ht analysis performed above should be redone with modihed values of R. 
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While this is beyond the scope of the present paper, it is left as a future task. 

The form ~ {xg)"^ for the 2-gluon density also neglects possible longitudinal correlations 
which should suppress the 2-gluon density at sufficiently large values of x. Also the possible 
difference of the fractional momenta of the fusing gluons is neglected. It can be argued, based 
on a simple picture of a Lorentz contracted proton and wavelengths of partons given by their 
inverse momentum, that gluons with x < l/{2mpR) ~ 0.1 overlap with any other gluons 
and thus cause finite nonlinear corrections. As the 1-gluon densities are already decreasing 
quite rapidly at a; >0.1, and as their \ogQ^ slopes become anyway small there, we have not 
attempted to build in any longitudinal correlations of the fusing gluons. In a larger system, in 
a big nucleus, this would be necessary, and sensitivity to such details should be studied. 

Another obvious improvement for the present analysis is to consider the GLRMQ terms 
added on top of the NLO DGLAP evolution. In that case, the zero (GTEQ6M) or negative 
(MRST2001) gluon distributions at small x may become larger than zero. 

In the present study we have neglected a possible but small contribution from the higher¬ 
dimensional gluon distribution Ght introduced in [8]. The distribution Get can be thought as a 
moment of /cT-dependent gluon distribution, and its upper limit estimated as Ght{x, Q^) ~ 
{k^)g{x, Q^) < Q^g{x, Q^). Therefore, this term should be less than the normal DGLAP term, 
and remain negligible. Further studies on this interesting question will follow [13]. 

Regarding the size of the nonlinear terms, it is quite interesting to notice that in our results 
at Qq = 1.4 GeV^, the GLRMQ terms in Eq. (1) make about 48 % of the full DGLAP-l-GLRMQ 
logQ^ slope of the gluon distribution at a; = 10“^, and still some 16% at a; = 0.01. The extent 
of nonlinearity decreases, as expected, with increasing Q^: at = 10 GeV^ the GLRMQ 
contribution to the total dxg/dlogQ^ is 26 % at a; = 10“^ and below 4% at a; = 0.01. At 
the lowest values of and x probed by the HI data [4], we clearly are at the borderline 
of the applicability of the approach, i.e. close to the gluon saturation region, where the next 
terms in the nonlinear evolution equations are becoming important [15, 16, 17, 18]. The further 
nonlinear correction terms obviously enter with an alternating sign. In the region where the 
GLRMQ term in Eqs. (1) and (3) becomes as important as the DGLAP term, inclusion of 
the further corrections should thus decrease the net correction. Therefore, the results of the 
current paper in the region of smallest x and studied can be regarded as an upper limit of 
the small-a; gluon distributions. 

In conclusion, we have studied the effects of adding the nonlinear GLRMQ corrections to 
the LO DGLAP evolution equations, and especially the interplay between the initial conditions 
and the nonlinearities. We use the PDF sets GTEQ5L and GTEQ6L as a baseline, and the 
recent DIS data from HI [4] as a constraint. We have shown that the agreement between the 
measured and computed structure function F 2 {x, Q^) can be improved at small values of x and 
while still maintaining the good fit to the data obtained in the global analyses at larger 
values of x and Q^. The nonlinearities slow down the scale evolution, so in order to recover the 
measured logQ^ slopes of F 2 {x,Q^) of the data, a larger small-x gluon distribution than that 
in GTEQ6L is needed. For the gluon distribution the nonlinear effects are found to play an 
increasingly important role at a; < 10“^ and < 10 GeV^. The nonlinearities, however, vanish 
rapidly at larger values of x and Q^. Gonsequently, contrary to GTEQ6L, the obtained gluon 
distribution at = 1.4 GeV^ shows a power-like growth at small values of x. Relative to the 
GTEQ6L gluons, an enhancement up to a factor 6 at a; = 10 = 1.4 GeV^ reduces to a 
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negligible difference at Q'^ > 10 GeV^. 
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Figure 1: The scale evolution of the structure function F 2 {x,Q‘^) of the free proton for fixed values 
of X (with constants added to separate the curves). The dashed curves show the LO DGLAP result 
from CTEQ5L [6], and the solid curve the result after the DGLAP+GLRMQ evolution when initial 
conditions taken from CTEQ5L at Qq = 5 GeV^. The data is from HI [4] and the error bars are 
statistical. 




2 


N(x) 



Figure 2: The goodness parameter of the fits of the computed T 2 (x, to the HI data, divided by 
the number of data points, as a function of the x of the data. The cumulative number of the data points 
is increasing to the left as indicated at the top of the plot: N{x = 0.2) = 2 and N{x = 3.2-10“®) = 133 
(see also Table 1). The curves are the LO DGLAP results from CTEQ5L (long dashed thick line) 
and CTEQ6L (dotted-dashed thick line), the DGLAP+GLRMQ result with the initial conditions at 
= 5 GeV^ taken from CTEQ5L (densely dotted) and from CTEQ6L (sparsely dotted), and, our 
set 1 (solid), set 2a (double dashed) and set 2b (short dashed) . 
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F2(x,Q )+c(x) 



3: As Fig. 1 but for LO DGLAP result from CTEQ6L (dotti 
+GLRMQ results with our set 1 (solid), set 2a (double dashed) and 






Figure 4: The parton distribution functions at = 1.4 GeV^ as obtained in the DGLAP anal¬ 
yses CTEQ5L [6] (dashed), CTEQ6L [3] (dotted-dashed) and in the present work based on the 
DGLAP-I-GLRMQ evolution. In our set 1 (solid), there is a finite charm contribution, while in the 
sets 2 (double dashed) charm is zero at this scale. Notice the enhancement from the GTEQ6L glue. 
The gluon distributions of set 1 and sets 2 are identical. 
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Figure 5: The dependence of the gluon distribution function at fixed values of x, from CTEQ5L [6] 
(dashed), CTEQ6L [3] (dotted-dashed) and set 1 of the present work (solid). Notice the logarithmic 
scales and the absolute normalization of the curves. 





